On the convergence of Filon quadrature

AbstractWe analyze the convergence behavior of Filon-type quadrature rules by making explicit the dependence on both k, the parameter that controls the oscillatory behavior of the integrand, and n, the number of function evaluations. We provide explicit conditions on the domain of analyticity of the integrand to ensure convergence for n→∞

FEM–BEM coupling for the large-body limit in micromagnetics

AbstractWe present and analyze a coupled finite element–boundary element method for a model in stationary micromagnetics. The finite element part is based on mixed conforming elements. For two- and three-dimensional settings, we show well-posedness of the discrete problem and present an a priori error analysis for the case of lowest order elements

Multiscale Petrov-Galerkin method for high-frequency heterogeneous Helmholtz equations

This paper presents a multiscale Petrov-Galerkin finite element method for time-harmonic acoustic scattering problems with heterogeneous coefficients in the high-frequency regime. We show that the method is pollution free also in the case of heterogeneous media provided that the stability bound of the continuous problem grows at most polynomially with the wave number k. By generalizing classical estimates of Melenk (Ph.D. Thesis 1995) and Hetmaniuk (Commun. Math. Sci. 5, 2007) for homogeneous medium, we show that this assumption of polynomially wave number growth holds true for a particular class of smooth heterogeneous material coefficients. Further, we present numerical examples to verify our stability estimates and implement an example in the wider class of discontinuous coefficients to show computational applicability beyond our limited class of coefficients

Variational Multiscale Stabilization and the Exponential Decay of Fine-scale Correctors

This paper addresses the variational multiscale stabilization of standard finite element methods for linear partial differential equations that exhibit multiscale features. The stabilization is of Petrov-Galerkin type with a standard finite element trial space and a problem-dependent test space based on pre-computed fine-scale correctors. The exponential decay of these correctors and their localisation to local cell problems is rigorously justified. The stabilization eliminates scale-dependent pre-asymptotic effects as they appear for standard finite element discretizations of highly oscillatory problems, e.g., the poor $L^2$ approximation in homogenization problems or the pollution effect in high-frequency acoustic scattering

Multiscale Partition of Unity

We introduce a new Partition of Unity Method for the numerical homogenization of elliptic partial differential equations with arbitrarily rough coefficients. We do not restrict to a particular ansatz space or the existence of a finite element mesh. The method modifies a given partition of unity such that optimal convergence is achieved independent of oscillation or discontinuities of the diffusion coefficient. The modification is based on an orthogonal decomposition of the solution space while preserving the partition of unity property. This precomputation involves the solution of independent problems on local subdomains of selectable size. We deduce quantitative error estimates for the method that account for the chosen amount of localization. Numerical experiments illustrate the high approximation properties even for 'cheap' parameter choices.Comment: Proceedings for Seventh International Workshop on Meshfree Methods for Partial Differential Equations, 18 pages, 3 figure

An adaptive variable order quadrature strategy

In this article we propose a new adaptive numerical quadrature procedure which includes both local subdivision of the integration domain, as well as local variation of the number of quadrature points employed on each subinterval. In this way we aim to account for local smoothness properties of the function to be integrated as effectively as possible, and thereby achieve highly accurate results in a very efficient manner. Indeed, this idea originates from so-called hp-version finite element methods which are known to deliver high-order convergence rates, even for nonsmooth functions

On thin plate spline interpolation

We present a simple, PDE-based proof of the result [M. Johnson, 2001] that the error estimates of [J. Duchon, 1978] for thin plate spline interpolation can be improved by $h^{1/2}$. We illustrate that ${\mathcal H}$-matrix techniques can successfully be employed to solve very large thin plate spline interpolation problem

Heat flow and calculus on metric measure spaces with Ricci curvature bounded below - the compact case

We provide a quick overview of various calculus tools and of the main results concerning the heat flow on compact metric measure spaces, with applications to spaces with lower Ricci curvature bounds. Topics include the Hopf-Lax semigroup and the Hamilton-Jacobi equation in metric spaces, a new approach to differentiation and to the theory of Sobolev spaces over metric measure spaces, the equivalence of the L^2-gradient flow of a suitably defined "Dirichlet energy" and the Wasserstein gradient flow of the relative entropy functional, a metric version of Brenier's Theorem, and a new (stronger) definition of Ricci curvature bound from below for metric measure spaces. This new notion is stable w.r.t. measured Gromov-Hausdorff convergence and it is strictly connected with the linearity of the heat flow.Comment: To the memory of Enrico Magenes, whose exemplar life, research and teaching shaped generations of mathematician

Finite Element Convergence for the Joule Heating Problem with Mixed Boundary Conditions

We prove strong convergence of conforming finite element approximations to the stationary Joule heating problem with mixed boundary conditions on Lipschitz domains in three spatial dimensions. We show optimal global regularity estimates on creased domains and prove a priori and a posteriori bounds for shape regular meshes.Comment: Keywords: Joule heating problem, thermistors, a posteriori error analysis, a priori error analysis, finite element metho

Flux norm approach to finite dimensional homogenization approximations with non-separated scales and high contrast

We consider divergence-form scalar elliptic equations and vectorial equations for elasticity with rough ($L^\infty(\Omega)$, $\Omega \subset \R^d$) coefficients $a(x)$ that, in particular, model media with non-separated scales and high contrast in material properties. We define the flux norm as the $L^2$ norm of the potential part of the fluxes of solutions, which is equivalent to the usual $H^1$-norm. We show that in the flux norm, the error associated with approximating, in a properly defined finite-dimensional space, the set of solutions of the aforementioned PDEs with rough coefficients is equal to the error associated with approximating the set of solutions of the same type of PDEs with smooth coefficients in a standard space (e.g., piecewise polynomial). We refer to this property as the {\it transfer property}. A simple application of this property is the construction of finite dimensional approximation spaces with errors independent of the regularity and contrast of the coefficients and with optimal and explicit convergence rates. This transfer property also provides an alternative to the global harmonic change of coordinates for the homogenization of elliptic operators that can be extended to elasticity equations. The proofs of these homogenization results are based on a new class of elliptic inequalities which play the same role in our approach as the div-curl lemma in classical homogenization.Comment: Accepted for publication in Archives for Rational Mechanics and Analysi